Direct measurement of the Q-function in a lossy cavity 
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Abstract 

We show that the Q-function corresponding to an electromagnetic field in a lossy cavity can be 
directly measured by means of a simple scheme, therefore allowing the knowledge of the state of 
the field despite dissipation. 
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The reconstruction of a quantum state is a central topic in quantum optics and related 
fields [f, 2]. During the past years, several techniques have been developed, for instance 
the direct sampling of the density matrix of a signal mode in multiport optical homodyne 
tomography [3], tomographic reconstruction by unbalanced homodyning [4], reconstruction 
via photocounting [5], cascaded homodyning [6] to cite some. There have also been proposals 
to measure electromagnetic fields inside cavities [7, 8] and vibrational states in ion traps [7, 
9] . In fact the full reconstruction of nonclassical states of the electromagnetic field [10] and of 
(motional) states of an ion [11] have been experimentally accomplished. The quantum state 
reconstruction in cavities is usually achieved through a finite set of selective measurements of 
atomic states [7] that make it possible to construct quasiprobability distribution functions 
such as the Wigner function, that constitute an alternative representation of a quantum 
state of the field. 

Nevertheless, in real experiments, the presence of noise and dissipation has normally 
destructive effects. Schemes that treat a lossy cavity have been proposed [8] that involve 
a physical process that allows the storage of information about the quantum coherences of 
the initial state in the diagonal elements of the density matrix of a transformed state. Here, 
we would like to re-take the problem of reconstruction of the cavity field as studied in [7] , 
but allowing the cavity to have losses (i.e. a real cavity). We then will show, that although 
it is not possible to reconstruct the Wigner function, it is still possible to recover whole 
information about the initial state via the Q-function. 

Let us first consider the ideal case of no dissipation. We consider then the Hamiltonian 
for the interaction between a quantized field and a two-level atom reads (we have set h = 1) 



H = ua^a H — ^-o z + A(a^<5"_ + aa + ), (1) 

where cv and a are the creation and annihilation operators for the field mode, respectively, 
obeying [a,^] = 1. <r + = \e)(g\ and <t_ = |<?)(e| are the raising and lowering atomic 
operators, respectively, |e) being the excited state and \g) the ground state of the two-level 
atom. The atomic operators obey the commutation relation [<t + ,<t_] = cr z . uj is the field 
frequency, u eg the atomic frequency and A is the interaction constant. When we have the 
condition on the detuning, 5 = oj eg — uj, 
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y > Vn+1 (2) 



2 



for any "relevant" photon number, we can obtain an effective interaction Hamiltonian in the 
dispersive limit (see for instance [12]) 



Hf = x tfaa z , (3) 

with x — A 2 /5. 

Let us consider the atom initially in the following superposition 

l^(0)> = -^(|e> + b» (4) 

and the state of the field to be arbitrary, denoted by \iPf(0)). If before the interaction we 
displace the unknown field by a quantity a, we have the state for the total initial wave 
function given by 

m^0)) = ^=(\e) + \g))\M^0)), (5) 

where \ip F (a;0)) = D(a)\tp F (0)) with D(a) = exp(aaJ — a*a) is the displacement operator 
with amplitude a. By calculating the average value of the electric dipole a x = a_ + <r + , we 
obtain [7, 13] 

oo 

(a x ) = °) cos(2 X mt) (6) 

m=0 

where the photon distribution P m (a;0) = \{^F(a; 0)|m)| 2 with \m) a number state. By 
choosing t — ir/ (2%), we end up with 

oo 

(* x ) = £^n(a;0)(-ir (7) 

m=0 

that, except for a factor of tc is the Wigner function [7, 14]. 

The above treatment does not consider dissipation, and it is the aim of this manuscript 
to study the case when we have a dissipative cavity. The problem of the interaction of a 
two-level atom with a quantized field in the dispersive regime in a cavity with losses was 
treated exactly by Peixoto and Nemes [12]. Here we will use a superoperator technique to 
solve it in an alternative way. 

In the interaction picture, and in the dispersive approximation, the master equation that 
governs the dynamics of a two-level atom coupled with an electromagnetic field in a high-Q 
cavity is 



±p = -i[H? f ,p]+£p 

where 

Cp = 2^apa) — '-fa^ap — '-fpa^a, 

and p the density matrix of the system. 
We define the superoperators 

Lp = —Ta^ap — pT^a^a, 

and 

Jp = l^apa} 

where we have defined 

f = liA + ix&z, 
with 1a = |e)(e| + \g)(g\- It is not difficult to show that 

[J,L]p = -S r Jp, 
where the superoperator §r is defined as 

S rP = tp + pf f . 

The solution to Eq. (8) subject to the initial state (5) is then given 

p{t) = e^ +J ^p(a; 0) = e Lt e^ J p(a; 0) 

where 

1 _ e -Srt 
f(t) P = — g P- 

and p{a; 0) = \ip{a; 0))(ip{a; 0)|. 



We need to operate the density matrix with the exponential of superoperators given 
above. It is not obvious how e^^ J will apply on (5), and therefore we give an expression for 
it 

e f(t)J m a )=jr ^&( a )\M0))(M0)\DK<*)] x (it) 

n=0 

This is, because / is an atomic superoperator, it will operate only on atomic states, and J 
will operate only on field states. It is not difficult to show that 

V n \^A)(M = [ (e + e) J x 1 A + 1 (2 e $)w J \e)(g\ + 1 (2 6 e)n j \g)(e\, (18) 
with £ = 7 + From (18), (17) and (15) we calculate (<7 X ) and obtain 

OO I /yi £ I OO 

<*«> = 2 E ^ E e ' 2 ^K fc l^( a )l^(o))l 2 + cc. 

m=0 ' ' ' fc=0 

, OO ry± £ OO / ,s . 

= V^E e " 2fc?tl!!i 7T^ : ( A; + Hp(0;a)|A; + m)+c.c. (19) 

2 z — ' m! * — ' fc! 

m=0 fc=0 

By changing the summation index in the second sum of the above equation, with M = m+k, 
we obtain 



<»■> = \ £ E ^Tj&^WOi «)i M > + < 20 > 

Tri=0 A// — m ^ ' 



m=0 ' ' M= 



Finally, we can start the second sum of (20) from M = (as we would only add zeros to the 
sum, because the factorial of a negative integer is infinite), and exchange the double sum in 
it, to sum first over m, which gives 



(a x ) = \ jr + Y 26t ) M (M\p(0;*)\M) + c.c. 

M=0 \ ? / 



(21) 



By defining 



and 



= /_, + e -^[sin(2x)-,co S (2r)n 

V V + e VT [cos(2r) + rj sin(2r)] / v ; 



■ V + e" 4 ^ + 2r)e- 2r > T sin(2r) \ 5 
^-1 J (23) 
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with t = yt and 77 = 7/x, we can have a final expression for (<j x ) 

00 

(^> = ^/i M cos(M^)(M|p(0;a)|M). (24) 

M=0 

Equation (24) in general differs from an s-parametrized quasiprobability distribution (see 
for instance [14]) because for 77 ^ 0, i.e. the dissipative case, we will have 9 ^ n and therefore 
cos(M9) 7^ e M , where e is a parameter. However, if we consider the case when /i — 0, the 
only term that survives in (24) is M = 0, and we obtain the Q-function. 

= (0|p(0;a)|0) = (a|p(0)|a) = 7rQ(a). (25) 

To look for the value that makes // equal to zero, we re-write, (23) for r = 37r/4 in the 
form 

3ir \r] - e~ 3vw/2 \ 
MT = T } = (1+^)1/2 > (26) 

and therefore, for 77 = exp(— 3rj7i/2) we obtain // = 0. The value for 77 may be easily obtained 
numerically with the result r/ M=0 ~ 0.274457. In Fig. 1 it is shown the behaviour of \x for 
different values of 77, for rj = (dashed line) // = 1 as expected and for 77 .274457 (solid 
line) it is seen that \x goes to zero at r = 37r/4. 

By direct measuring a quasiprobability distribution, namely the Wigner function in the 
lossless case [7], i.e. for only one parameter (7 = 0), or the Q-function in our case for a range 
of parameters (7 .274457%), one can obatin complete information about the state of the 
field as quasiprobability distributions contain complete information of the density matrix 
[15]. Moreover the Wigner and Q functions are related by integral or differential equations, 
for instance [14] 

W(a) = e~s&&Q(a). (27) 

In conclusion, we have solved the dispersive interaction between a quantized elecromag- 
netic field and a two level atom in the case of a real cavity (subject to losses) by utilizing 
superoperator techniques. We then have used those results to show that even in the dis- 
sipative case we can still obtain information about the initial cavity field by means of the 
Q-function (unlike the non dissipative case [7] where it is reconstructed the Wigner function). 
Both functions, being quasiprobability distributions contain complete information about the 
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state of the cavity field. In this way, we have been able to extend the range of parameters in 
which complete information may be obtained in CQED, from 7 = to 7 =rs .274457%, i.e. 
by doing a right tuning complete information of the cavity field may be obtained despite of 
cavity losses. 

FIG. 1: We plot \i as a function of r for r\ = (dashed line), 77 = 0.025 (dashed-dotted line) and 
rj « .274457 (solid line) 
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